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Color superconductivity in high density QCD exhibits the color-flavor locked (CFL) phase. To ex- 
plore zero modes in the CFL phase in the presence of a non-Abelian vortex with an SU(2) symmetry 
in the vortex core, we apply the index theorem to the Bogoliubov-de Gennes (BdG) Hamiltonian. 
From the calculation of the topological index, we find that triplet, doublet and singlet sectors of 
SU(2) have certain number of chiral Majorana zero modes in the limit of vanishing chemical poten- 
tial. We also solve the BdG equation by the use of the series expansion to show that the number of 
zero modes and their chirality match the result of the index theorem. From particle-hole symmetry 
of the BdG Hamiltonian, we conclude that if and only if the index of a given sector is odd, one 
zero mode survives generically for a finite chemical potential. We argue that this result should hold 
nonperturbatively even in the high density limit. 

PACS numbers: 21.65.Qr, 11.27.-fd, 12.38.-t, 25.75.Nq 
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I. INTRODUCTION 

Zero modes around a vortex in a superconductor have 
been attracting much current interest, providing us with 
intriguing notion of the Majorana fermion. Vortices in 
a. Px -\- ipy superconductor yield zero-energy Majorana 
bound states ]H which obey non-Abehan statistics 0- 
II]. Genericahy, these degenerate states are expected to 
be quite useful in constructing a fault-tolerant quantum 
computer 0, [lo| , and the idea for topological quantum 
computation utilizing non-Abelian anyons on a topolog- 
ical state of matter has been developed [ll|. Majorana 
states have been predicted in many other systems such as 
a surface state of a topological insulator with the prox- 
imity effect of an s-wave superconductor [l^ or of a fer- 
romagnetic insulator [Tsj . an s-wavc superfluid of ultra- 
cold atoms [l^, the superfluid '^He B-phase [l5|, etc [l6l - 
[23} . Not only in condensed matter physics but also in 
high density QCD with color superconductivity [24l[25|. 
Majorana zero modes along an Abelian vortex [26| and 
along a non-Abelian vortex [23| in the CFL phase have 
been discussed, and new non-Abelian statistics has been 
derived [2^. Fermionic zero modes in the presence of 
topological background have also been explored in the 
standard model of the electro- weak interaction [29l - l33 |. 
Recently, a classification scheme for zero modes associ- 
ated with topological defects has been proposed [H, [s^] , 
which is a generalization of the top olog ical classification 
for insulators and superconductors [STllssj. 

Various methods of counting zero modes in such sys- 
tems have been proposed [^, [39l - l48l. Among them, appli- 
cation [131 of the index theorem j39l l40j is unique in that 
it does not resort to any approximations. The index the- 
orem has a long history and has elucidated, for example, 
topological characteristics of anomalies in gauge theories 



[49| . The index theorem applied to the present problem is 
a variant, dealing with a fermion coupled to a condensate 
with a nontrivial topological defect. This rigorous theo- 
rem claims equivalence between following two quantities: 
One is associated with the number of zero modes which 
reflects the analytical property of a differential equation. 
This can be obtained if one indeed solves the eigenvalue 
equation [s^, [Hi • The other is associated with a topolog- 
ical invariant of the order parameter with a defect. The 
index theorem thus relates the zero modes with the topo- 
logical configuration of a defect. It has been applied to 
a fermion describing a topological superconductor in 3D 
[5^ , and generalized to a Zq index theorem for a fermion 
without chiral symmetry [53]. 

In this paper, we apply the index theorem to the CFL 
phase of a color superconductor with a non-Abelian vor- 
tex A = diag(AQ, Aq, Ag) with an SU(2) symmetry in 
the vortex core, where Aq and Ag have winding q and 
Q, respectively. We find that triplet, doublet and singlet 
sectors of SU(2) have q, Q, and 2Q fermionic zero modes, 
respectively, propagating along a vortex line in the limit 
of vanishing chemical potential. These modes can be re- 
garded as one-dimensional chiral Majorana fermions due 
to particle-hole symmetry. We obtain these results from 
the calculation of the topological index as well as from 
the analysis of the normalizability of the zero mode wave 
functions, and find complete agreement between them. 
We argue that when the chemical potential is switched 
on, no (or more precisely, even) zero modes are expected 
in the case of an even index, whereas one (or more pre- 
cisely, odd) zero mode survives in the case of an odd 
index. This holds in each sector of SU(2), and is due to 
spectral symmetry between positive- and negative-energy 
states of the BdG Hamiltonian. Therefore, it should be 
extrapolated even into very high density limit nonpertur- 
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batively. 

This paper is organized as follows. In the next section 
im we derive effeetive Hamiltonians for triplet, doublet, 
and singlet states of SU(2) from the whole BdG Hamil- 
tonian. In Sec. IIII[ we give a brief review of the index 
theorem relevant to the present problem, and calculate 
the topological index for all sectors of SU(2). In Sec. IIV| 
we switch to solving the BdG equation, and count the 
zero modes and their chiralities directly, from which we 
know the analytical index. We show that these two in- 
dices indeed coincide. In the final section |Vl summary 
and discussions are given. In AppendixEl we show that a 
singlet zero mode found in [27| together with a new solu- 
tion diverges in a short distance and is non-normalizable 
although it was shown to converge in a large distance. In 
Appendix[Bl we give a brief review of the way to solve the 
differential equations of zero modes in a series expansion. 



II. EFFECTIVE HAMILTONIAN 

Under a non-Abelian vortex with an SU(2) symmetry, 
all eigenstates in the CFL phase are decomposed into 
multiplets of SU(2). Therefore, for later convenience, we 
reduce BdG Hamiltonian into smaller ones which act on 
each multiplet. 



A. BdG equation and vortices 

We start with a generic Hamiltonian with a given pair- 
ing gap of superconductivity. 



H = 



(fa 



(2.1) 



where a, /? and i, j are color and flavor indices, respec- 
tively, C = i'j'^j'^ is the charge conjugation matrix, and 



A^^ is a generic gap function with A^^ 
Fermi statistics. In the last line, we have introduced the 
Nambu representation 



,Q/3 



A^*" due to 



(2.2) 



where {ipfY = ^-0"^ = i^'^ipf* is the charge conjugation 
of '(/'f , and the total Hamiltonian density Jitot is defined 
by 



(2.3) 
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In what follows, the hat means the Nambu representa- 
tion with particle-hole symmetry, where and (V'f )'^ 
are referred to as a particle and a hole, respectively. This 
Hamiltonian density leads us to the BdG eigenvalue equa- 
tion. 



(2.5) 



Counting the zero modes of this equation is the main 
interest of the present paper. Here, the zero mode means 
a massless mode propagating along a vortex which will 
be specified momentarily. 

Now we take three flavors (i,j =) u,d, and s, and 
assume the following gap function in the CFL phase 



a: 



a/3 



(2.6) 



which belongs to 3 representation with respect to SU(3)c 
as well as SU(3)f- Under the action of {e^^ ,Uc,Uf) £ 



U(l)i 
as 



SU(3)q X SU(3)p, the gap function transforms 



A^ ^ e" 



The gap function takes a form of 



(2.7) 



(2.8) 



in the ground state, where Acfl is a constant. In 
this state, the so-called color-flavor locking symmetry 
SU(3)c+F, which acts on the gap as 



(C/AC/t) 



(2.9) 



with U € SU(3), is preserved. 

We consider a vortex state in the CFL phase in this pa- 
per. The minimal winding non-Abelian vortex, denoted 
by All in [131 , is given by 



diag(Ao,Ao,Ai)^ 



(2.10) 



where Ag(r, I 



|Ag(r)|e"''' stands for a gap function 



with a g-vorticity which obeys the boundary conditions 
|Ao(c»)| = |Ai(oo)| = Acfl and |Ao(0)|' = |Ai(0)| = 0. 
This vortex carries a 1/3 quantized circulation of U{1)b 
and a color magnetic flux is confined inside its core 
which is neglected in the BdG equation. The CFL 
symmetry SU(3)c+F in Eq. (j2.9p is spontaneously bro- 
ken into SU(2)c+F X U(1)c+f around the core of the 
vortex. Consequently. Nambu-Goldstone bosonic zero 
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modes CP2 ~ SU(3)c+F/[SUr2lc+F X U(1)c+f] appear 
inside the core of the vortex [63 jSSl - ISSt . Properties of 
this vortex has been extensively studied in the Ginzburg- 
Landau model [53 - l59| , and in the BdG equation [2g . 
On the other hand, a non-Abelian vortex, denoted by M2 

Ai; = diag(Ai,Ai,Ao), (2.11) 

with the same boundary conditions on the diagonal el- 
ements Aq as in Mi, preserves the same symmetry 
SU(2)c+F X U(l)c+F in the vortex core This vortex 
carries a 2/3 quantized circulation of U{1)b- 

More generally, in this paper, wc consider a bound 
state of q Mi and Q M2 non-Abelian vortices, for which 



the gap function takes the form of 

A!;==diag(AQ,AQ,A,). (2.12) 

The case of g = 1 and Q = is the non-Abclian Mi 
vortex above, while q = and Q = 1 corresponds to the 
non-Abelian M2 vortex. When q Q, this vortex breaks 
SU(3)c-HF to SU(2)c+fx U(1)c-hf- Therefore, the BdG 
eigenstates in p.Sp are classified as the SU(2) multiplcts. 
To derive an effective Hamiltonian in each multiplct is 
the main task in this section. 

For a generic gap function in Eq. (|2.12p , wc adopt the 
basis V' ~ {ur, dg, s;,, dr, Ug, Sr, Ub, Sg,db)'^ to write down 
the BdG equation in Eq. such that [gO, HH 
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(2.13) 



where Aq is defined by 



A - ( ^9757" 



(2.14) 



and Aq similarly. 

Let us now derive an effective Hamiltonian which acts 
on each multiplct of SU(2)c+f- To this end, it is conve- 
nient to write the particle states as a 3 x 3 matrix via 

= (V').a, 



Ur Ug Ub 

dr dg db 

Sr Sn Sb 



(2.15) 



where (a = 1, • • • , 8) is the generator of SU(3) normal- 
ized by liToTb = ^a6/2 defined through the Gell-Mann 
matrices as Ta — Xal'2- and Tg = (l/V^ • 1). This equa- 
tion implies that the particle states of the CFL phase 
can be regarded, if they are in the ground state specified 
by Eq. (|2.8p . as the adjoint representation of SU(3)c-i-f, 
Y^a=i ^oTa and an invariant state of a U(l) subalgebra 
^fgTg. The non-Abelian vortex in Eq. (|2.12p tells that 
we can reduce this representation into those which are 
irreducible under SU(2)c+f generated by the subalgebra 
(Ti , T2 , Ta ) . Under such subalgebra, (Ti , , ) them- 
selves transform as the adjoint (spin-1) representation. 



(T4 + iT5,Te + iTr) and (Tg - iTr, T4 - iTg) transform as 
the spin- 1/2 representation, and Tg and Tg transform as 
spin-0 representation. Therefore, the states (|2.15p are de- 
composed into one triplet, two doublets, and two singlets 
such that 113 



*t = Ug{Ti + m) + [Ur - dg)T3 + dr{Ti - 1T2), 

*di = Ub{Ti + iTs) + 4(T6 + iTj), 
*d2 = Sg(Tfi - iTj) + Sr{Ti - iTs), 

Ur+dg - 2sb 
^sl — 7= J8, 



x/3 



V6 



(2.16) 



Likewise, the hole states can be written as 



(2.17) 



since they are complex conjugate of the particle states, 
ip'^ = i^'^tp* . Therefore, corresponding SU(2)c+f multi- 
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plets are 



= + iT^) + K - ^9)^3 + - zT2), 
= d^(r4 + in) + ulin + ^tt), 



V3 



n 



*S2 = 2- 



V6 



-n. 



(2.18) 



It is now easy to see that the BdG equation in (|2.13p is 
decoupled into smah pieces acting on the above muhi- 
plets, as win be done below. Among them, our primary 
concern is the singlet states. 



B. SU(2) doublet sector 

The two doublet states, whose spaces are spanned by 
Mf), dfc, Sg, and Sr, are coupled only to the gap Aq. Equa- 
tion (|2.13p shows that these are decoupled into {sr,Ub) 
and {sg,db)- Considering Hq in Eq. (|2.4p and Aq in 
Eq. (|2.14p . we see that among the states {sr,Ub), for ex- 
ample, {sr,ul) and (ub,s'^) are decoupled, and (sg,dh) 
likewise. This means that the Hamiltonian for {sr,Ub) 
and (sg,dfc) is block-diagonal, and BdG equation is de- 
composed into 



nj 

no 



E 



(2.19) 



where 



((p, rj^) ^ (s,,, <), {ub, <), (sg, dg), (4, s^). (2.20) 



It should be noted that the Hamiltonian in Eq. (|2.19p 
can be simply given by 



Ho 



(2.21) 



Thus, this can be regarded as the effective Hamiltonian 
in the doublet states. We here mention that from Eqs. 
(Pl^ and (I23H1), the SU(2) quantum number n reads 
n = -1/2, 1/2, 1/2, and -1/2, respectively, for the 
paired states {sr,ul), (Mb,s^), (sg,(ig), and (dh,Sg). This 
implies that the initial particle and hole paired states, 
e.g., (sr,s5;), switch their partners to (sr,uj;) and so on 
in order to be the simultaneous eigenstates of n St,s well as 
the Hamiltonian. In this representation, the new paired 
states are related by 



(2.22) 



where anti-unitary particle-hole transformation C is de- 
fined by 



^ " ' i^^ 



i"f 



K, = 1, 



(2.23) 



with complex conjugation K. Correspondingly, the effec- 
tive Hamiltonian (|2.2ip has indeed particle-hole symme- 
try, 



C(Ho-Aq)C-i = -(?{o-Aq), 



(2.24) 



as it should do. Although the BdG equation in the dou- 
blet states is described by single equation (|2.19p . the 
above particle-hole symmetry implies that if ((/?, 77^^) is 
regarded as particle states with an energy E, correspond- 
ing hole states are given by (77, ip"^) with energy —E. The 
Hamiltonian (|2.2ip . which appears in Eq. (|2.19p . can 
thus be regarded as an effective Hamiltonian involved in 
the SU(2) doublets with correct SU(2) quantum numbers 
{J2a=i 2 3^a'^3) ^'^'^ with corrcct particle-hole symme- 
try described by C. 

This Hamiltonian can be further decomposed into 
right- and left-handed modes, i.e., the eigenstates of 75 
with eigenvalues 75 = ±1. Let {tp,ri'^) be a paired SU(2) 
doublet state defined in Eq. (|2.20p . and let 95 = ((/?+, (/?_) 
and 77^^ = ('7l,'7+) be the decomposition into right and 
left modes, where -|- and — mean the right and left modes, 
respectively. We use the following 7 matrices in this pa- 
per, 



-a \ I I 

75 



-1 



(2.25) 

Then, the BdG equation (|2.19p can be decomposed into 



-i(T • V - M -Aq W V?+ 



-A* icT■^ + ^l \ It 



= E 



(2.26) 



and 



icT -XI — [I A(; 



in -S/ + fj, J \ 



(2.27) 

These equations define effective Hamiltonians ?ij which 
act on SU(2) doublets with 75 = ±1 chirality, 



Ho- Aq 



-Ho A, 



Tit 



(2.28) 



where 



A± = ± 



-i(T ■ V 



ia- ■ V 



-1 



(2.29) 



When we calculate the index in the next section, it is 
convenient to define new T matrices such that 

P=crJ®CT-\ (j = 1,2,3) 



= 1 (g) 



(2.30) 
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which satisfy the relation {r'',r^} = 2(5'^^. Then, the 
Hamihonians in the above can be written as 



where 



-A*n icr ■ V 



(2.31) 



with j = 1,2,3, a = 4, 5, and 

Cj)= (04,05) 

= (-ReAQ,ImAQ) 



(2.32) 



(2.33) 



The Hamiltonian in Eq. ([QT]) is the final expression 

of the efi'ective Hamiltonian in the doublet states, and Tid 
in Eq. ()2.32p plays a crucial role in the index theorem in 
the next section. 

This Hamiltonian belongs to the gene ric universality 
class of superconductors, class D (s^Hlll, since it trans- 
forms 



(2.34) 



under a reduced particle-hole transformation C defined 
by 



1, 



(2.35) 



which is associated with Eq. (|2.23p , acting on the spaces 
of right or left modes separately. 



C. SU(2) triplet sector 

The triplet states couple only to the gap Aq, and when 
q = 1, the explicit solutions of the zero modes have al- 
ready been obtained in |27[ . Equations (|2.16p and (|2.18p 
tell that three states Ug, dr, and ii. 
in this sector. As to the T3 
follows from Eq. (|2.13p that they obey the same BdG 
equation as in the doublet sector if the gap function is 
replaced by A^, and therefore, an efi'ective Hamiltonian 
for paired state (dr,Wg) and {ug,d'^) is given by 



dg are involved 
±1 states, Ug and dr, it 



"Ho — A„, 



(2.36) 



in the same sense as in Eq. (|2.2ip . As to the T3 = state 
iir — d, 

Eq. (|2.13p . the following equation 



g, it is easy to derive, from the BdG equation in 



{Ho - Aq){Ur ~ dg) ^ E{Ur - dg). 



(2.37) 



whose Hamiltonian is just the same as the one in 
Eq. (|2.36p . Therefore, the Hamiltonian (|2.36p can be re- 
garded as an effective Hamiltonian in the triplet sector. 
Chiral decomposition is also carried out in the same way 



as in the previous subsection, and the final Hamiltonian, 
which is referred to as "H^, is given by Eqs. (|2.3ip and 
(|2.32p with a different order parameter 

= (04, 05) 

= (-RcA,,ImA,). (2.38) 
This Hamiltonian also belongs to class D. 



D. SU(2) singlet sector 

While the doublet and triplet sectors have rather sim- 
ple Hamiltonian basically equivalent to single flavor case 
[23, the singlet sectors couple to gaps with different 
winding, which could yield nontrivial zero modes. Actu- 
ally, the two singlet states are coupled together through 
Aq and Aq, whose BdG equations are derived from Eq. 
such that 



(Ho + Ag){Ur + dg) + 2AqS6 = E{Ur + dg), 

HoSb + Aqiur 



dg) = Esb- 



(2.39) 



Here, note that any linear combination of singlet states 
is always singlet, and hence Ur -\- dg and Sf, are two sin- 
glet states. Let {{ur +dg)/^/2, Sb)"^ be an extended wave 
function. Then, we can define a Hamiltonian which leads 
to Eq. (fOg)) such that 



Uo + Aq V2Aq 
V2Aq -Ho 



(2.40) 



This Hamiltonian can be regarded as an effective Hamil- 
tonian in the singlet sector. Chiral decomposition similar 
to Eq. (|2.28p leads to the Hamiltonian 



(2.41) 



where and Ag (and A^ similarly) are defined by 
Eq. (|2.29p . Let us define new F-matrices 




= K 
F^ = 1( 
F^ = 1( 



a 



8) 1, 
)1, 



(j = 1,2,3) 



r =l'g)a^ (g) a^, 
which satisfy the relation {F^, F"} = 26^"- 

plp2p3 = jl^cr^^l, 



(2.42) 
Note that 



p4p5p7 
p5p6p7 ^ , 
p6p4p7 ^ ^ 



' a 
-il (g) 



1, 



(2.43) 
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Then, by the use of the new T matrices, the above Hamil- 
tonian can be written as 

Ht = ±Hs + ifiT^T^T^, (2.44) 

where 

/ -ict W Aq ^/2AQ \ 
y ^ a; jrr-V V2A*Q 
V2Aq -i(T-V 
V V2A*Q ia W J 

= -iV^dj + r"(^a + ^i^^^'^VaF'TT^. (2.45) 

Here, two kinds of indices take j = 1,2,3 and a = 4, 5, 6, 
stands for the completely antisymmetric tensor with 
^456 = +1^ and 

(f> = (04,05,06) 

= (V2ReAQ,iReA„-ilmA,), 
ip = (■04, V'5, V'e) 

= (y2Im Aq, ^Im Ag, ^Re A,). (2.46) 

The Hamiltonian thus obtained in (|2.44p acts on the 
two singlet states with 75 = ±1, respectively. 

This Hamiltonian also belongs to the same universality 
as the previous multiplets, class D, since it has particle- 
hole symmetry 

CntC-^ = -nt, (2.47) 
with a reduced particle-hole transformation 

C = iT^r^T'^K, & = 1. (2.48) 

III. THE TOPOLOGICAL INDEX 

In the previous section, we have decomposed the to- 
tal Hamiltonian into several pieces acting on irreducible 
multiplets under SU(2)c+f with a definite chirality of 
75 = ±1. These Hamiltonians, denoted as H.^j with 
m —A^ t, and s, belong to universality class D of topo- 
logical superconductors. We will omit the subscript m 
for simplicity in this section below. As discussed by Teo 
and Kane [Sa, [s^ , a line defect in class D allows in gen- 
eral protected zero modes of Majorana type propagating 
along a defect. 

A. Zero modes along a vortex 

To explore such Majorana zero modes, we assume that 
the vortex has a cylindrical profile described by the cylin- 
drical coordinates (r, 6*, z) such that 

A,(r,0,z) = |A,(r)|e'«''W, (3.1) 



where generic angle function Qq{0) has a winding q, 
0g(27r) = Og(0) -|- 27rg. We also assume for the radial 
part that |Ag(c>D)| ~ Acfl = const. > as well as 
|A,(0)| = if (7 7^ {\Aq{0)\' = for q = 0). To in- 
vestigate the number of gapless modes, we first consider 
the case of the zero chemical potential, fi = 0. For any 
multiplets, the Hamiltonian is then given by ±1-1 , as can 
be seen from Eqs. (|2.3ip and (|2.44p . Since the ± sign 
has nothing to do with the zero modes, only the Hamil- 
tonian H will be henceforth investigated. The cylindrical 
symmetry in Eq. p.ip enables us to separate the motion 
to the z-direction from others, 

H{r,9,k,) ^n^{r,e) +k^T\ (3.2) 

where the Fourier transformation has been made for the 
z-direction. Suppose that the eigenvalues for T-L^ are 
obtained such that 'H±^e = E^e- Then, by the use of 

{?ii,r3} = 0, (p^O), (3.3) 

we see that nj_(T^'^E) = -£;(f^^£;). We will refer to 
this property as F'^-chiral symmetry of TL^l . It thus turns 
out that the eigenfunctions of H arc linear combinations 
of ^E and F'^^I'b, and hence, their eigenvalues are those 
of 

f E \ 
\k, -E J' 

namely, iy/k'^ + E'^. It follows that an i? = state 
of 'H± yields a zero mode of H: Let 'I'o be a zero- 
energy state of H±. Due to F'^-chiral symmetry, we can 
choose it as a simultaneous eigenstatc of F''. Suppose 
F^^Poi = ±*o±- Then, we see that ?i*o± = ifc^^oi- 
Here, note that the phase of the wave function 5'o±(^, 
can be chosen such that 

C^o±ir,0) = i'o±ir,0), (3.4) 

because of particle-hole symmetry in Eqs. (|2.34p or 
(|2.47p . This tells that the present mode along a vortex 
line is a chiral Majorana zero mode. Therefore, the num- 
ber of such Majorana modes of H, denoted by NuiTi), 
is the same as the number of zero-energy bound states of 
'H±, denoted by No{'H±), which are classified by integers, 
i.e., Z. Namely, 

Nuin) = No{ni_), iii = 0). (3.5) 

On the other hand, if a nonzero chemical potential fi is 
switched on, it breaks the F'^-chiral symmetry, since the 
chemical potential term in Eqs. (|2.3ip and (|2.44p . which 
will be referred to as T^^i, is not anti-commutative with 
F^. Let us write the Hamiltonian as 

n^{r,9,k,) = ±n^{r,9)+n^+k,T^ 

^ V ' 

= n^{r,e)+k,T\ (3.6) 
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and let us consider the zero-energy states of ii^{r,9). 
When fj, is small, we can treat the term "H^ as a pertur- 
bation to the unperturbed degenerate zero-energy states 
oi'H± [15 . We then expect that without any special sym- 
metries, even number of zero modes of 'H± become gener- 
ically nonzero-energy modes with energies ie^-, whereas 
odd number of zero modes allow at least one unpaired 
zero mode. It may sometimes happen by chance that 
several zero modes occur if we fine-tune parameters of 
the model, but the evenness or oddness of this number 
should be invariant. We can see this alternatively from 
the symmetry property of the model: Particle-hole sym- 
metry guarantees that the spectrum is symmetric with 
respect to the zero energy. This implies that if we con- 
sider the spectral flow as a function of /x, the number of 
zero modes modulo 2 should be invariant even towards a 
very large values of /i, 

NoiUf) ^ No{H^) mod 2, {^i^0). (3.7) 

Therefore, the above number is a topological invariant 
characterizing the class D superconductors. 

Next let us consider the zero modes of "H*. As has 
been discussed, when /i = 0, we have simultaneous eigen- 
states of T-L^ and T^, implying they are also eigenstates 

of for any kz- This is due to F'^-chiral symmetry 
when /.J = 0. Contrary to this case, a zero-energy state 
of is not an eigenstates of any longer when /-t 7^ 0, 
which makes it difficult to find out an eigenstate of 
for generic kz- However, at least when kz = 0, it is 
indeed an eigenstate of il^ with the zero energy, and 
therefore, we can claim the existence of a gapless mode, 
although we cannot know its exact dispersion relation. 
In Ref. [13], an effective theory for this gapless mode 
has been derived. Even with a finite chemical potential, 
the model has particle-hole symmetry and we can choose 
wave function with Majorana-like condition. We thus 
conclude that the number of gapless Majorana modes of 
for /i ^ is equivalent to the number of the zero- 
energy states of for /i = modulo 2, 

Nuin^) ^ Noin±) mod 2, (m^O). (3.8) 



In this sense, even if we are interested in CFL phase 
realized at very high density QCD, it is important to 
study the case of the zero chemical potential. Therefore, 
we will concentrate on the number of zero-energy states 
of 'H_L in each multiplet. 

Now we switch to the discussion on the index of the 
Hamiltonian. When /i = 0, we can define the following 
index, 

mdH±^ N+{-H^_) - N^iHi_), (^ = 0), (3.9) 

where N±{'H±) is the number of zero-energy states of 
with the definite F'^-chirality, F"^ = ±1, respectively. 
This is actually possible, since Eq. p.3p guarantees that 
zero-energy states can be simultaneous eigenstates of F'^. 
Since No{'H±) = N+{n±) + N^{n±), we have No{n±) = 
ind {T-L±) mod 2. Together with (|3.7|) . we finally reach 

A^mCH=^) = ind CH^l) mod 2, {fi ^ 0). (3.10) 

Although wc have defined the number of zero modes mod- 
ulo 2, the typical number is and 1 when the index is 
even and odd, respectively. We expect that any other 
numbers can appear only by chance or due to some sym- 
metries. 



B. Index Theorem 

The index theorem states that the index defined above 
in Eq. (|3.9p can be written by the topological invariant. 
To see this, let us rewrite the index as follows [iol. |4^: 

indH± = lim TrF^^^ , (3.11) 

where Tr stands for the trace over the 2D coordinate 
space as well as over the F-matrices. In the derivation of 
the index theorem, a central role is played by the axial 
vector current defined by 



J^x,m,M) = lim trF^F^ . ]S^(x~y), 

= lim trF^F^' (m±) I ^r-^ — | S'^(x - y), (3.12) 

where we have introduced a Pauli-Villars regulator with a large mass parameter M to make the current well-defined. 
After all calculations, we should take the limit M — >■ 00. Suppose that H± = —iT^dj + (57^ j_, as in the case above. 
Then, the divergence of the current yields the index such that 

d,J^{x,m,M) = 2 lim trF^ - — 5'{x^y). (3.13) 

y^^ \n^, +7772 -^2 _j_j^p / 
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Therefore, we reach 



1 



ind'H± = c+ hm ^ (h e,jJ\x,m,M)dx^, (3.14) 



where 



c = hm Tr F 



(3.15) 



It is known that c becomes the Chern number associated 
with the gauge potential [i^. In the present case, it 
vanishes, since we have neglected a color magnetic flux 
confined in the vortex. 



C. Doublet and triplet states 

Let us first review the calculations of the topological 
index, i.e., the right- hand-side of Eq. p.l4p in the simpler 
cases of the doublets and triplet, which is basically the 
same as the model already studied [13, 13] . 

In both cases of doublets and triplet, the effec- 
tive Hamiltonian in Eq. p.l4p is constructed by 
Eqs. (|232l) and with Eq. (|233| for the doublets 

and with Eq. (|2.38|) for the triplet. The momentum rep- 
resentation for Eq. p.l2p yields 



_fk_ 



1 



ikx 



^2 



(3.16) 



where j = 1,2, and contributions from the regulator have 
been neglected, since the above current is well-defined, as 
seen below. This is due to the fact that for j = 1,2 and 
the order parameter V^fj^a are ant i- commutative. How- 
ever, it should be noted that if one considers more generic 
systems, it plays an important role generically. For con- 
venience, especially for comparison with the singlet sec- 
tor, we introduce the following notations. 



where 



= (?ii +m2)e^'^^ = X - A, (3.17) 



A = iT^r''dj<l>a + 2ikjdj + a| (3.18) 

with k"^ = Y^j=i ^1 a^*i 'P'^ = Sa=i 4''a- Here (pa is given 
in Eq. ([05]) for the doublet sectors and Eq. (P?^ for 
the triplet. Expansion {K - K)-^ = E„(-^"^A)"iC-i = 
X;„ A'7is:"+i leads us to 



(fk 1 

(2^^ 



YtrF'T^(i-H± + tr*fc^)A". 



(3.19) 



This expansion is quite useful, since the index theo- 
rem p.l4p needs the current only at |a;| = r — > 00: 



The assumption \4>\ Aqfl = const, at r — J- 00 
tells that dj(j)a ^ 0{r^^). Therefore, we see that the 
terms n > 2 in Eq. p.l9p vanish when substituted into 
p.l4[) . On the other hand, the 71 = term also van- 
ishes by the trace. Therefore, only 77 = 1 term can con- 
tribute to the index, and we finally reach, by the use of 
trF3FJF''F''F'' = {2i)^ei^e'''' , 



dPk 1 



(27r)2 ^2 
(27)2e^-^e«''0,9,0, 
1 



-trF^r^irVb^r^r^Sf^fc + o{r-'^) 



(Pk 1 



(27r)2 



(3.20) 



At the infinity r 00, we see <p/\4>\ = (— cos 8q, sin 6q) 
for the doublets, and 4>/\4>\ = (— cos 9^, sin 8^) for the 
triplets. Taking the limit 777 0, we finally have 



This result claims that the 2D Hamiltonian ii± allows at 
least Q (g) zero modes for each doublet (triplet) state, 
and therefore, the 3D Hamiltonian with a vortex line 
fi in Eq. p.2p has at least Q {q) gapless modes along 
the vortex. When the chemical potential is taken into 
account, it follows from the discussion in Sec. IIII Al that 
the doublet (triplet) allow a zero mode if Q {q) is odd. 
We have found that the non-Abelian Mi vortex {q = 1 
and Q = 0) has one triplet zero mode as explicitly found 
in [231, while the non-Abelian Mi vortex (g = and 
(5=1) turns out to have one doublet zero mode. 



D. Singlet states 

So far we have calculated the index of the 2D Hamil- 
tonian 7i_L for both the doublet and the triplet states, 
which are basically the same Hamiltonian already stud- 
ied. On the other hand, the Hamiltonian "H^ for the 
singlet states, which is defined by Eqs. (|2.45p and p.2p . 
is rather complicated and quite unique, since two kinds 
of different vortices, Ag and Ag. are involved simultane- 
ously. In the limiting case where Aq — 0, the Hamilto- 
nian is decoupled into upper and lower halves. The 
upper is completely the same as the triplet Hamiltonian, 
whereas the lower is just a free fermion Hamiltonian. The 
index of the former and the latter is q and 0, respec- 
tively. However, once the nontrivial Ag is switched on, 
it is quite interesting to ask which winding, A^ or Ag, 
controls the zero modes of the Hamiltonian. 

Although the calculations are parallel to IIII C[ the ex- 
istence of the terms composed of three F-matrices in 
Eq. (j2.45p makes them quite complicated. First, note 
that 
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= K + T-A, 

with j = 1,2 and a,b,c ~ 4, 5, 6, where (j)a and ■i/'a arc those in Eq. p.46p for the singlets, and 



(3.22) 



(3.23) 



with (IP = 4't ^i^d V'^ = X]a=i '/'a- Then, the expansion of the kernel {i-L\ + vn?) ^ is also quite useful. 



(2^ 



e-*'=^trr'^r(m_L) 



1 



Akx 



j (^*'^'^' + '^^'^ + " (i/2)i?'"'=V^ar^rT7] 



A' + r - A 



E 

ri=0 



<fk 

(2^ 



trr^r* [r^(aj + ik^) + ir^^^ - (i/2)£;°^'='(/>ar''rT^] 



A - r 



A 



K -T 



(3.24) 



Here, we have also neglected the regulator with mass M . In the last equation, it is enough to calculate only the n — 1 
term for the same reason as in Sec. IIII CI After taking the trace, we have 



<Pk 



(27r)2 (i4r2 -r2)2 



{ [K^ - 2K{(b^ + V^2) + T^]^.1f,^ + 2i^(02 _ V'')^(0 ■ ^)} , (3.25) 



where (t>-ip — Sa=i 4>aipa and (j) ■ ''di tp = (j) ■ diip — {dgcp) ■ ip- By the use of r2 = 4 (^cp'^ip'^ — {(j) ■ ip)^) and the momentum 
integrations 



d^k 



K 



1 



1 



(27r)2 (A:2 - r2)2 StT (02 + V'2 + m2)2 - 4(l/'2?/'2 - (0 • V)2) ' 

dPk K'^+T'^ _l 02 + V^+m2 

(27r)2 {K^ - r2)2 ^ 4^ (02 + ^2 ^2)2 „ 4(^2^2 _ (0 . ^)2) ' 

I 



we end up with 



TT (02 + -02 + 77^2)2 _ 4(^2^2 _ (0 . ^)2) ' 

(3.26) 



When the limit m ^ is taken, the following two cases 
(1. Aq ^ and 2. Aq = 0) should be considered sepa- 
rately. 



1. Generic case: Aq 7^ 

In this generic case, including the case of Q = (but 
Aq=o 7^ 0), we can take the limit m directly in 



Eq. (13^. Then, we have 



, ,,2 (02-^2)g^(^.^) 

(02 _ ^2)2 _(_ 4(^ . ^)2 



7r|$|2 



where 



V;2,20.V') 



1*1 ^(02 _,/,2)2_,_ 4(^.^)2 

= (cos26q, sin26Q). 



(3.27) 



(3.28) 



This equation implies that the index of the Hamiltonian 
is determined solely by the vortex Aq. Moreover, the 
index is twice the winding number of Aq: 

indTii = 2^ idede{2QQ) = 2Q. (3.29) 
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This is the central result in the former part of the present 
paper: The topological index of the singlet states is 2Q 
in generic cases. Therefore, in the typical case of g = 
1 and (5 = (but Aq=o 0), we conclude that the 
index is zero, implying no zero modes generically. Of 
course, we cannot deny the existence of even number of 
zero modes, since the index is not the number of zero 
modes. However, in such cases, those zero modes are 
unstable, and small perturbations yield gaps for them 
[4^ . In this sense, we expect in general no zero modes 
in the singlet states. Actually, in the next Sec. IIV[ we 
solve the differential equation for zero modes directly, 
and reach the same conclusion. 

We have thus concluded that the non-Abelian Mi vor- 
tex (9 = 1 and (3 = 0) and M2 vortex ((7 = and Q = 1) 
has no singlet zero modes, but it seems to be inconsis- 
tent with the result in [23| , in which an asymptotic form 
of a singlet zero mode at large distance from the vor- 
tex core was given for the Mi vortex. In Appendix |X] 
we will show that the given asymptotic solution actu- 
ally diverges at short distance r — >■ around the vortex 
core. We will present another solution which is also well- 
defined at r — ^ 00. However, is is diverges as well at short 
distance. Therefore those are non-normalizable modes, 
which the index does not count. We thus conclude that 
no contradiction exists. 



2. Exceptional case: Aq = 

Next, let us consider the exceptional case, Ag = 0. 
Although we do not need any concrete calculations to 
know the index as discussed already, we show this case 
below to check the validity of our calculations. When 
Aq = 0, a naive limit m ^ makes the denominator 
of Eq. p.26p vanish. However, note that t/)^ = -0^ and 
(f> ■ ip = Q when Aq = 0. Thus, we can take the limit 
m 0, 



.1^ 



2 rv?<^ ■ di -tp 



7r|(/)|2 



where 4>/\(j)\ — (0, cos 9g, — sin 6g). Therefore, 



(3.30) 



(3.31) 



We have thus obtained the expected index: As discussed, 
when Aq = 0, the Hamiltonian decouples into two pieces. 
One is the same as the Hamiltonian for the triplet states, 
and the other is a free fermion Hamiltonian. It is obvious 
without any calculations that the former and the latter 
give the index q and 0, respectively. Once a nonzero Aq 
is introduced, however small it may be, the index jumps 
from q to 2Q, as seen above. 



IV. ANALYTICAL INDEX - 
ZERO MODES 



COUNTING 



So far we have calculated the topological index for the 
doublet, triplet, and singlet states as a winding number 
of the pairing gap function. In this section, we explore 
normalizable solutions of the equation 



(4.1) 



To this end, we assume Qq{0) = q9 and 8q(6') = Q9 
in this section. The merit of such an analysis is that it 
informs us of the number of zero-energy states and their 
F'^-chiralities, and separately. 



A. Doublet and triplet states 

For doublet and triplet states, the Hamiltonian H± is 
defined by Eqs. (|2.32p and p.2p . In a suitable basis in 
which = diag(l, 1, -1, -1), Eq. (|i?T|) becomes 





n±.+ 



0, 



where ± denotes the r-^-chirality, = ±1, and 



—id± 
|A,(r)|- 



|A,(r)|e'9» 

idzr 



(4.2) 



(4.3) 



with d± = di±id2 = e^'^(9r±^9e). This is the equation 
for the triplet, and if q is replaced by (5, it becomes for 
the doublets. We set the (r, 6')-dependence of the wave 
functions such that 



m± 



a™(r)e±™» 



(4.4) 



where to is a quantum number associated with the angu- 
lar momentum. Then, the equation for the radial part is 
given by 



ar r 



(4.5) 



where i' = (a™,/?™)' 



M± = diag (to, —to — 1 ± g), 
|A„ 



|A„ 







(4.6) 



Counting zero-energy states for this model has been car- 
ried out by Jackiw and Rossi [Slj . As they have assumed, 
it may be physically natural to consider the asymptotic 



gap function at r 
generically 



as |A„ 



'-I'l. We here assume 



n = ^f7„r", 



(4.7) 



n=0 



11 



where f2o = 0, when g ^ 0. From a mathematical point 
of view, the difference between \Aq\ ~ r'*' and \Aq\ ^ r 
docs not affect the number of zero-energy states. How- 
ever, the latter is given by log-corrections in the series 
expansion of the solutions. Sec Appendix |B] for details. 
Since Eq. (|4.5|) is composed of two first-order differential 
equations, general solution has two parameters: One of 
these parameters is determined by the normalization of 
the wave function. As a result, the general solution has 
one free parameter. Let tp{r) be the general solution of 
Eq. (|4.5p . We analyze the asymptotic behavior of ^^{r) 
below. 

First, let us start with the consideration on the behav- 
ior of tplr) at r — > oo, where we can neglect the term 
1/r in Eq. (|4.5p . The gap function is constant there, 
l^gC'')! ~^ ^CFL (> 0). Then, the set of equations 
(|4.5p allows two independent solutions ^ e^^"^ , where 
A± = zbAcFL are eigenvalues of fi. It follows that the 
general solution ili{r) becomes at r oo a linear com- 
bination of such exponentially decreasing and increasing 
functions ^^{r) — X]i=± '^^^'^''^i where Oi is a constant 
vector. For the wave function to be normalizable, we im- 
pose one condition on the wave function that g+'^cFL'- 
should vanish, a+ = 0. This means that one free pa- 
rameter mentioned above is determined completely. One 
normalizable solution of Eq. (|4.5p thus corresponds to 
one zero-energy solution for T-L^. 

Next, let us consider the normalizability of the wave 
function at r 0. For in Eq. (|4Jl) . Eq. (|4?5l) can 



be solved by the series expansion techniques. Details are 
given in Appendix [B] It turns out that two particular so- 
lutions are allowed in this limit as well, which are written 
in power series of r such that 

^(') = V^Wr"' + V^^Vi^"'^' + • • ■ ' (4-8) 

for i = 1,2, where the leading power Ui is one of the 
diagonal elements of M, i.e., rii = m and 712 = — m — l + g 
for F"^ = +1 and ni = m and 712 = —m — 1 — q for F'^ ~ 
— 1. It should be noted that the general solution ip{r) 
becomes a linear combination of (r) in the limit r — > 
and that ip has no free parameter any longer. It follows 
that both ■f/''"^'' and -0'^' should be normalizable. The 
normalizability of ■0''^ in the limit r — >■ is guaranteed 
by Hi > 0. Namely, < < 9 ~ 1 for F'^ = -f 1 and < 
771 < —q — 1 for F"^ = —1. Accordingly, (1) when I < q, 
N+ = q and 7V_ = 0, (2) when q ^ 0, N+ = = 0, 
and (3) when g < — 1, iV_|- = and iV_ ~ ~q. Thus, we 
end up with indH^ = A^-i- — N_ = q'va any cases above. 
We have thus reproduced the index in Eq. p.2ip indeed. 



B. Singlet states 

For singlet states, the Hamiltonian "H^ is defined by 
Eqs. (|2.45p and p.2p . In a suitable basis, the zero mode 
equation becomes Eq. (|4.2p with 



-id^ 





V V2|AQ(r)|e- 



|A,(r)|e^9e 

V2|AQ(r)|e"^ 






^/2|AQ(r)|e-^0« 
~id± 




V2|AQ(r)|e'«» 





idm 



(4.9) 



We set the (r, 
that 



dependence of the wave functions such 



/a™(7^)e±™^ 



(4.10) 



Define ip = (a,„, /?,„, 7,„, (5,„)"^. Then, the equation for 
the radial part is given by the same as Eq. (|4.5p . where 



AI± = diag (m, —m — 1 ± q, m ± (Q — q) , —m — 1 ± Q), 
/ \A,\ V2\Aq\\ 

■ (4.11) 





\V2\Ac 



^/2|Ac 
V2\Aq\ 




We also assume the series expansion Eq. (|4.7p for fi. 
Counting the zero-energy states will be carried out simi- 



larly to the previous case. From Eq. (|4.1ip . it is obvious 
that when |AQ(r)| = 0, Eq. (|4.5p is reduced to two pieces; 
One corresponds to the triplet states studied above, and 
the other is just free fermions. Therefore, in the following 
analysis, we assume |AQ(r)| 7^ generically. 

Equation (|4.5p is now composed of four differential 
equations. Therefore, in the present case, the general 
solution includes four parameters, and one of them is de- 
termined by the normalization of the wave function. The 
general solution has thus three free parameters. This is 
the difference between the present singlet case and pre- 
vious doublet or triplet cases. Let ip{r) be the general 
solution of Eq. (|4?5|) . 

Let us first consider the asymptotic behavior of ip{r) at 
r — >• 00, where the 1/r term can be neglected in Eq. (|4.5p . 
We assume that two gap function approach the same 
constant, |Aq(c!o)| = |Ag(oo)| = Aqfl = const. (> 0). 
Then, ipi^r) becomes a linear combination of the four in- 
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dependent solutions ip — ^ ■ a^e'^''' with A±i = ±Acfl 
and X±2 = ±2AcFL, where a, are constant vectors. For 
the wave function to be normahzable, we should impose 
two conditions a+i = and a-|_2 = 0. Therefore, the 
general solution has still one free parameter. This is in 
sharp contrast to the previous case in which the general 
solution has no free parameters at this stage. 

Let us next consider the behavior of 'ijj{r) at r — > 0. 
In this region, we can solve Eq. (|4.5p by using the power 
series expansion of the wave functions. For details, see 
Appendix|B] We show there that the set of four equations 
allow four special solutions (|4.8p whose leading power 
{i = 1, • • • ,4) is given by the diagonal elements of M± 
defined by Eq. (|4.1ip . It follows that, near r 0, we 
have 



(4.12) 



4=1 



where c; is a constant. When we consider the normal- 
izability of the wave function ip{r), one free parameter 
plays a crucial role, since one of Ui can be negative. In 
what follows, we restrict our discussions to F"^ = +1 case, 
for simplicity. 

Suppose Q > 0. Then, the four powers are given by 
ni = TO, n2 = —TO — I + q, = m + Q — q and ~ 
—TO — 1 + Q for F'^ = +1 states. We should have some 
cases separately. 

(I) Case ofQ <q. 

(1) If all Hi are non-negative, < ni,2,3,4, we 
have q — Q<m<Q — 1. In this case, all 
are normahzable at r — > 0. 
This implies that we do not need to impose 
any conditions of normalizability on ip{r). Re- 
member here that the general solution ip{r) 
has one free parameter. Therefore, ipir) is a 
linear combination of two independent solu- 
tions. In other words, this case allows degen- 
erate two solutions for each quantum number 
TO. Therefore, we have 



2{2Q 




q) (g<2Q) 
(2Q < q) 



(4.13) 



where, the factor 2 above means two indepen- 
dent solutions. 

(2) If ni < 0, < ^2,3.4, no to is possible. 

(3) If 712 < 0, < ^1^3. 4, no TO is possible. 

(4) If < 0, < '^1,2,4, we have 



< TO < 



q-Q-l {q<2Q) 
0-1 (2Q < g) 



(4.14) 



Since 713 < 0, this particular solution should 
be eliminated from tp{r) at r — )■ given in 



Eq. (|4.12p by imposing the additional condi- 
tion C3 = 0. This determines one free param- 
eter which 'ip{r) has, and normahzable ip{r) 
thus obtained gives one solution for each quan- 
tum number to. Therefore, 



q-Q iq<2Q) 
Q {2Q < q) 



(5) If 7i4 < 0, < 711^2,3, we have 
<m <q 1 < 5) 



(4.15) 



(4.16) 



For the same reason as in the case of (4), for 
each allowed to, one zero-energy state is given, 
and the number of the zero-energy states are 
thus 



q-Q {q<2Q) 
Q {2Q < q) 



(4.17) 



Summing up iV_|_ over the cases from (1) to (5), we 
conclude that the number of zero-energy state is 
given by iV_|- = 2Q in the case oi Q < q. 

(II) Case ofq<Q. 

(1) If < ni,2,3,4, we have < m < q — 1. Since 
no additional conditions of the normalizability 
need to be imposed, there are two degenerate 
solutions for each to. Therefore, we have 



2q {0<q) 
iq<0) ■ 



(4.18) 



Note here the factor 2, as discussed in (I) (1) 
(2) If Til < 0, < 712,3,4, we have 



q — Q < m < 



-1 (0 < q) 
q-1 {q<0) 



(4.19) 



Since each allowed to gives one zero-energy 
state, and the number of the zero-energy state 
is given by 



N. = 



Q-q iO<q) 
Q {q<0) - 

(3) If 722 < 0, < Til, 3, 4, we have 



<TO<0-1 ^^-'^^ 

-""-^ ^ iq<0) 



(4.20) 



(4.21) 



The number of the zero-energy state is 



Q-q iO<q) 
Q {q < 0) 



(4.22) 



(4) If 723 < 0, < 711,2.4, no TO is possible. 
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(5) If < 0, < '^1,2.3, no TO is possible. 

It turns out that the number of zero-energy states 
is also given by = 2(5 in this case. 

So far we have counted the zero-energy states with — 
+1 in the case of Q > 0. In a similar way, we immediately 
see that there are no solutions allowed in the = — 1 
states. Thus we conclude that iV+ = 2Q and = 
when Q > 0. On the other hand, in the case of Q < 0, 
we can show that iV+ = and -/V_ = 2\Q\. These results 
match precisely the topological index given by Eq. (|3.29p . 

In the above analysis, we have used just the leading 
power of ipi (.T^) ■ In Appendix [B] wc discuss the higher 
orders of the series expansion and show how to construct 
the solutions. This formulation may be interesting on 
its own right, in particular, in that we need to introduce 
(Inr)" (n < 3) terms to obtain four independent solu- 
tions. 



V. SUMMARY AND DISCUSSIONS 

Non-Abelian vortices exist in the CFL phase of color 
superconductivity. We have calculated the topological 
and analytical indices for fermion zero modes of non- 
Abelian vortices A = diag(AQ, Ag, A^) which preserve 
SU(2)c+F symmetry in their core, where Ag and Aq 
have winding Q and q, respectively. Because of SU(2)c+f 
symmetry unbroken in the core of the vortex, the BdG 
Hamiltonian is decomposed into some sectors according 
to the irreducible representations of SU(2)c+f; triplet, 
doublet and singlet sectors. In the limit of the zero chemi- 
cal potential /i = 0, the topological indices for the triplet, 
doublet and singlet zero modes of SU(2)c+f have been 
calculated to be q, Q and 2Q, respectively. We have also 
analyzed the BdG equations, and from the normalizabil- 
ity of the wave function we have obtained the same in- 
dices of the zero modes. Because of particle-hole symme- 
try, these zero modes can be regarded as chiral Majorana 
fermion modes. For non-zero chemical potential 7^ 0, 
we have expected that the triplet or doublet zero modes 
exist only when g or Q is odd, respectively, whereas no 
zero modes in the singlet states. Technically speaking, 
the derivation of the indices of the triplet and doublet 
states is essentially the same with that of one component 
(iol HI], but that of the singlet states is quite nontriv- 
ial; the index of the singlets is q if we set Ag to zero by 
hand, but it jumps to 2Q once non-zero Ag is turned on, 
however small it may be. 

In particular, as the most fundamental non-Abelian 
vortices, the non-Abelian Mi vortex (q = 1 and Q = 
0) carrying 1/3 quantized circulation of [/(1)b, has one 
triplet, no doublet and no singlet Majorana zero modes. 
Although the number of triplets agrees with the previous 
result in that of singlet does not seem so at first 
glance. However we have shown in Appendix |X] that the 
asymptotic zero mode at large distance from the vortex 
core found in [27| diverges in the core. The index docs 



not count such a non-normalizable mode, and hence, no 
inconsistency exists. On the other hand, the non-Abelian 
AI2 vortex (g = and Q = 1) carrying 2/3 quantized 
circulation of U{1)b, has no triplet, one doublet and no 
singlet Majorana zero modes. 

Before closing this paper, several discussions are ad- 
dressed here. 

We have calculated the indices for non-Abelian vor- 
tices of a composite of q Mi and Q M2 types, where 
the SU(2)c+F xU(l)c-i-F symmetry remains in the core 
of vortices. Generalization to the case with three com- 
ponents having all different winding numbers is straight- 
forward, in which case only U(1)q^p symmetry remains. 

Although we have considered (composite) vortices 
placed in the same position, they can be separated if they 
are composite. In fact, a non-Abelian M2 vortex can be 
decomposed into two non-Abelian Mi vortices winding 
in different components. It is an interesting problem to 
see how the indices of SU(2) multiplets change from M2 
to Mi's in this process. Here, the positions of the vor- 
tices play a crucial role as a new degree of freedom. One 
doublet zero mode of one M2 vortex will disappear, as it 
separates into two Mi vortices, and when they are well 
separated in position, a triplet zero mode will appear in 
each of two Mi vortices. They are at the zero energy 
when two Mi vortices are infinitely separated. Thus, 
the doublet zero modes are lifted to non-zero modes in a 
process of decomposition of one M2 vortex into two Mi 
vortices, whereas new zero modes appear in turn in the 
triplet states. It implies the existence of attractive and 
repulsive forces between two Mi vortices, mediated by 
doublet and triplet (and maybe singlet) fermion modes, 
respectively. 

A new non-Abelian statistics of non-Abelian Mi vor- 
tices with triplet zero modes was studied j28j . in which 
it was shown to be written as a tensor product of that of 
Abelian vortices in a chiral p-wavc superconductor [3| and 
the Coxtcr group. Our result implies that non-Abelian 
M2 vortices have doublet zero modes, which may give 
another new non-Abelian statistics. 

Non-Abelian vortices should exist in the core of neu- 
tron stars if the CFL phase is indeed realized there. They 
will constitute a vortex lattice because of the rapid rota- 
tion of the stars [H^, HH . It is known that fermion zero 
modes on vortex lines significantly change the transport 
properties of quasi-particles (63j . Fermion zero modes 
found in this paper will be essential in the study of dy- 
namics of neutron stars, which may hopefully gives an 
observational signal. 

Finally, the CFL phase admits a variety of topolog- 
ical defects: domain walls, Sk yrm ions |64| . (confined) 
monopoles [59| and instantons [65|. A generalization to 
the index theorems in the presence of these topological 
defects will be an interesting work to be explored. It also 
counts edge modes at the phase boundary [6l|, be- 
tween the CFL phase and the hadron phase, if it exists. 
Also, the inclusion of strange quark masses (H^j and an 
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extension to the other phases such as 2SC phase [25| will 
be important in application to more realistic situation. 



Acknowledgments 

This work was supported in part by Grants-in-Aid 
for Scientific Research from JSPS (Nos. 20340098 and 
21540378) for Fukui and Grants-in-Aid for Scientific Re- 
search (Nos. 20740141 and 23740198) from the Ministry 
of Education, Culture, Sports, Science and Technology- 
Japan for Nitta. This work was also supported by the 
"Topological Quantum Phenomena" (No. 23103502 for 
Fukui No. 23103515 for Nitta) Grant-in Aid for Scientific 
Research on Innovative Areas from the Ministry of Edu- 
cation, Culture, Sports, Science and Technology (MEXT) 
of Japan. Yasui is supported by the Grant-in- Aid for Sci- 
entific Research on Priority Areas "Elucidation of New 
Hadrons with a Variety of Flavors" (No. 21105006) from 
the MEXT of Japan. 



Appendix A: Singlet solutions in the asymptotic 
form at large r 

The asymptotic form of a singlet zero mode of a non- 
Abelian Mi vortex at large r was given in [l^. This is 
actually a candidate of zero-energy states, which is ex- 
ponentially decreasing function at r — > oo. However, the 
asymptotics at r — is also crucial for the normalizabil- 
ity of wave functions. To investigate the behavior of this 
solution extrapolated into a small r region, we need to 
solve the BdG equation numerically. 

For numerical calculations with high accuracy, it is de- 
sirable that we have two independent solutions at r — > oo. 
A linear combination of them is then a general solu- 
tion including a parameter. Starting with such a generic 
wave function, we can investigate precisely the behavior 
of the numerical wave function extrapolated into r — > 
by changing the parameter of the generic wave function 
at the boundary r — )- oo. 

To this end, we will present not only the one shown in 
Ref. [131 but also another asymptotic solution at r — > oo 
in this Appendix, and report the result of the numerical 
computation of the BdG equation. 

We show the explicit form of the wave functions of the 
singlet solution for the right mode (75 = +1) in the Weyl 
representation 



where 



and 



Sb- 



(Al) 



(A2) 




(A3) 



(A4) 



for 



= 1) and dg {i = 2), and 



.h{r)e~' 









(A5) 



(A6) 



forsft. Here we have the relations /i (r) — —g2ii"), gi{r) ~ 
fjir), Mr) = -gi{r), gi{r) = /2(r), Mr) = -52 (r) and 
g3{r) = fsir) from the Majorana condition. At large r, 
both |Ao| and |Ai| become a common constant |A| given 
in the bulk state, |Ao|, |Ai| |A| = Acfl(> 0) at 
r — > 00. Then, with an approximation of small |A| and 
large ^, we find asymptotic forms of the wave functions 
with a condition of the convergence at large r. The first 
solution is 



Mr) = AAe-l^l'-/2jo(/ir), 
g,{r) = AAe-l^l'-/Vi(/ir), 



(A7) 
(A8) 



{i = 1, 2) and 



Mr) 
Mr) 



-:^e-IAk/Vo(H, 



e-l^l'-/Vi(H, 



(A9) 
(AID) 



with a normalization constant M . This is the solution 
given in the previous work [27| . As a second solution, we 
find a new asymptotic solution, which was not considered 
in the previous work. 



Mr) = ^^'e-\^\^'^^-{^,r)\M^^r) 

(Ji(H^o(/ir) - M^lr)Nl{^xr)), (All) 
g[{r) = AA'e"l^l'-/2l ( _ ^rJo(MO + Ji(H), (A12) 
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{i = 1, 2) and 

-7r(H' (^2Jo(H + ^Ji(h) 
(Ji(H^^o(H - Jo(H7Vi(H)}, (A13) 
g^(r) = AA'e-l^l'-/2_^{ (2 + |A|r) Jo(/ir) 

(Ji(HA^o(H - Jo(H7Vi(H)}, (A14) 

with a normalization constant Af'. It should be empha- 
sized that these asymptotic solutions are correct only at 
large r, at which |Ao| and |Ai| are constant. However, 
these solutions may be divergent in small r in general, 
because |Ai| becomes zero at r = 0, and we find that 
it is the case. In order to see the behavior at small 
r, we have solved numerically the BdG equation with 
assuming an approximate r-dependence of Ao(r) and 
Ai(r) (for example, Ao(r) = Acfl — const. > and 
Ai(r) = AcFLtanh^r with a coherence length ^) with 
starting from large r in which these asymptotic solutions 
arc given as a boundary condition. As a numerical result. 



we have found that the wave functions become divergent 
at r = for any linear combinations of the first and the 
second asymptotic solutions as a boundary condition. It 
means that these asymptotic solutions are not normaliz- 
able in the whole range of r. Therefore it is confirmed 
that there is no normalizable zero-energy state in the sin- 
glet, consistent with the result of the index theorem in 
the text. 



Appendix B: Analytic solution of the differential 
equations 

We have counted the zero modes based on a simple 
assumption that the leading power of the particular so- 
lutions for Eq. (j4.5p near r '--^ is given by the diagonal 
elements of the matrix M± in Eqs. (|4.6p and (|4.1ip . In 
this Appendix we show it indeed possible to find such 
solutions. 

We assume that the differential equation has a power 
series solution 

oo 

i^ = r^^7^A+nr". (Bl) 

Putting this as well as the Taylor expansion of Q in 
Eq. (|4.7p into the differential equation Eq. (14. 5p . we can 
write the differential equation as 



A-l 



^ i (A + n + 1 - il/)VA+„+i + J2 ""-'^ 



.n+\ 



(B2) 



n=0 



/=0 



where M± has been denoted simply as M. This gives 
(A - M)yjx = 0, (B3) 

n 

{X + n + l^ M)tPx+n+i + J2 ^n-ii'x+i = 0, (B4) 

(=0 

with n > for the latter. From Eq. (|B3p . it follows that 
A should be one of the diagonal elements of M. For the 
largest value of A, Eq. (|B4[) can be solved recursively, 
whereas for smaller A, we could meet the difficulty, since 
at a certain n, det(A + n + 1 ~ M) = 0. 

1. Case of triplet and doublets 

Without loss of generality, we can set 

M — diag(p, s), p> s, (B5) 

where p and s are integers. As mentioned above, the 
leading power A is given by A p and X = s. Corre- 



spondingly, Eq. (jB3[) gives two solutions V'^*'' (* = 1:2) 
as 

If p = s, we see det(n + A + 1 - A/) 7^ for n > 0. We 

can find V'i+A-i-i recursively by (jB4[) . We thus obtain two 
solutions. Henceforth, we assume p > q. 

a. Generic case of \Aq\ ~ r 

First of all, let us consider the most generic case and 
give the solutions of Eq. (|4.5p just by assuming Eq. (|4.7p . 
For X = p, higher (n > 0) can be obtained recur- 

sively by Eq. (|B4p . and the solution is given by 

00 

n=0 
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For the smaller power A = s, the recursion relation is 
ill-defined at n=p — s— l>0, since det{p — M) = 0. 
This has a close relationship with the well-known fact 
that the Bessel function J„ and J_„ are not independent 
when n is an integer, and the Neumann function Yn in 
addition to the Bessel function J„ can be regarded as two 
independent solutions of Besscl's differential equation. In 
the present problem, the standard technique of solving 
differential equations leads to the following ansatz 



2. Case of singlet 

This case is rather complicated, since Eq. (|4.5|) is com- 
posed of four coupled equations. We first rearrange the 
differential equation Eq. (|4.5p for the singlet case to sat- 
isfy 



M — diag(rii, 712, 71.3, 71.4), 
ni> n2> ny,> 714, 



(B12) 



ijj^'^^r) = ^ i^n+s'''"^'^ + a,4'^^\'r) Inr, (B8) ^iid introduce a set of normalized eigenvectors ipnl {i 



where a is an unknown constant to be determined. Sub- 
stituting this into Eq. (j4.5p modifies Eq. (jB4[) such that 

n 

{s + n + l- M)0^+, + ^^n-i41i + "V-l+Ui = 0' 

1=0 

(B9) 

where V'i+s+i = for 77 = 0, • • • ,p — s — 2 is assumed. 
The difficulty at7i=p— s — 1 due to det(p — M) = is 
avoided by introducing a parameter a. We can determine 
both a and ipp by this equation. For details, especially 
for the uniqueness of the solution, sec the discussions in 
App. IB 2 bl Once these arc obtained, we solve V'i^^ 



+ S+1 



for n > p — s recursively. We thus obtain the second 
solution (|B8p . Although it contains Inr, the expected 
leading power behaviors near the origin is not altered. 



h. Case o/|A,| ~rl«l 

So far we have shown that if the diagonal elements of 
M are different integers, the solutions include generically 
a ln7' term. We will show below that if we assume |Ag| ~ 
r'''! , the normalizable solutions of Eq. (j4.5|) do not include 
such a term, and they can be given by purely power scries 
of r. 

We restrict our discussions to F^ = -1-1 and g > 1, for 
definiteness. Then, 



M ~ — diag(?77, q 



1), 



(BIO) 



and normalizability requires < m < g — 1. We explore 
the case m > q — m — 1 only. The solution for the larger 
A = 777 is given by (|B7p . For the smaller A = — 777 — 1 , 

we can choose ip'f^m-i to be the second one in Eq. (|B6I) . 



Eq. (|B4p gives recursively V'n+g-m for < 77 < 2777 — q. 
The problematic equation at n — 2m — q can be satisfied 

(2) 

simply by choosing ipm ~ 0, since couplings with fig only 
appear in (|B4p for n > q, whereas < 27T7 — q < q — 2. 
Eq. (jB4p then determines ipn+q~m for n > 2m — q, giving 
the second solution 



E,/,(2) „"+g-m-l 
V-'n+g-m-l' 

n=0 



(Bll) 



Therefore, in this case, there are no terms including Inr. 



I,... ,4)by 



(1) 



7/;(4) 







(B13) 



The solutions of Eq. (|4.5p near the origin change their 
forms if any of 77^ are equal. Such degeneracies are clas- 
sified into eight different cases. We consider these sep- 
arately. In particular, when there are more than two 
different diagonal elements, there appear (Inr)" correc- 
tions with 77 > 2, which is a new feature absent in the 
previous triplet and doublet cases. 



a. Til = 712 = 7I3 = 774 



We define p = 771. Eq. (|B3I) gives four independent 
solutions Vp = V'p^ (i = 1, • • • ,4) for A = p. Eq. ((B4)) 



then recursively determine ipp+n+i = V'p+n+i- We thus 
obtain four solutions given by 



+p' ' 



(i = l,---,4). (B14) 



n=0 



h. 7ll = 712 = 713 > 7I4 

Let us define p = 771 and s = 774. Then, three solutions 
{i = 1, 2,3) w ith A = p can be obtained as in the 
case of App. IB 2 al The fourth solution for A = s can be 
assumed to be 

00 

= ^ + ^ a^^-'^^^{r)\nr, (B15) 

n=0 1=1,2,3 

where a^^'*^ are some constants. The recursion relations 
for the coefficients can be found by inserting Eq. (jBlSp 
into Eq. They read 



(s - M)V'i'*) = 0, 

n 

(s + 77 + 1 - M)^itl+i + E 



(B16) 



s+l 



1=0 



E «('''Vlt+i=0, (" = 0,1,...) (B17) 



i=l,2,3 



where V'i+n'''' = for n = 0, 1, • • • ,p — s — 2 is as- 
sumed. Eq. (|B16p is automatically satisfied by the fourth 

vector in Eq. (|B13p . Eq. (|B17p determines 

V'p-i recursively since det(s + n + 1 — M) ^ for 
n = 0, 1, • • • ,p - s - 2. Eq. (|B17| for n = p - s - 1 
needs some care because of det(p — M) = 0. It gives 
Qi(4:2) without ambiguity, whereas ^p^'^ cannot be de- 
termined uniquely since we can freely modify V'p^^ by 
Si=i 2 3 , where Ci are arbitrary constants. It can 
be seen that such an arbitrariness does not give any new 
independent solutions and can be removed simply by as- 
summg ^PI*'^ oc i.e., the fourth in Eq. ((BT3)) . Once 
'0p^'' and a^"*'*' are obtained, we can find V'i+„+i for ^ P- 

We shall be brief since we can infer the solutions for 
any given M in a manner similar to these. 



c. ni = n2 > ns = ua 

We define p = rii and s = 71,3. Two solutions ■0^'^ 
(i = 1,2) with \ ~ p can be obtained as in the case of 
App. IB 2 al Other two solutions V'*-*' (« = 3, 4) with A = s 
can be written as 

00 

= ^ ?/'i^,r"+" + ^ a^^J'V'^^'^ Inr. (B18) 

n=0 j = l,2 

We can find all the coefficients ■i/'i+s > 0) and a'^*'-'^ 
as in the case of App. IB 2 bl 
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e. ni > 722 = W3 = n4 

This case is similar to App. iBlbl or lBTcl We can find 
solutions of the form 

00 

00 

= ^ V^^_|,r"+" + a^'-^V'^^^ Inr, (z = 2, 3, 4) 

n=0 

(B20) 

where p = ni and s = 77.2. 



/. ni > 712 = 713 > 714 

The solutions can be found in the form 

00 

00 

= ^ V;,(;|,r"+*^ -f a^^'i^V'^^ Inr, (i = 2,3) 

ri=0 

= ^ V4'|tr"+* + fa(4'i) Inr + ^^(^^Danr)^') 

n=o V 2 y 

^ V^*' In ^, (B21) 

i=2,3 

where p = 77,1, s = 712 and t = 714. 



Til = 712 > 713 > 714 



g. Til > 712 > 713 = 714 
The solutions can be written as 



Defining p ^ ni and s = 713, we obtain three solutions 
t/i^') (i = Ij 2, 3) as in the case of App. IB 2 cl The fourth 00 
solution takes the form ^(i) = ^ Vi+p^"'*'^, 



n=0 

00 

=EV'a-"-^*+ E (a^^''^lnr + i6(^-^)(ln.)^) V'^ ^''^ = E ^it^"^^ + -^"'^^^'^ l^^' 

n=0 i=l,2 ^ ^ " ° 

+ a(4.3)V;(3)lnr, (B19) = ^ ^r"+* + fa^'^^hnr + -b^'^^\\nr)A ^^^^ 

. , ,1, 1 , x9 -Ha^^'^V^^^lnr, (i = 3,4) (B22) 

where t = n^. It should be stressed that (mr) term 

appears. In cases where M has more than two different 

diagonal elements, there appear higher Inr corrections. where p = tt-i, s = 712 and t — n^. 
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h. ni > 712 > ns > n4 



where p = ni, s = 712, t = and k = 714. 



(1) ^n+p 



n=0 



n=0 



+ (^a(4^2) Inr + l6(4,2)(^^^)2^ ^(2) ^ „(4,3)^(3) j^^^ 



(B23) 



In any cases, (Inr)" terms appear in nonleading order 
of r. We thus conclude that the solutions behave as 



^W(^) ^ (* = 1,--- ,4). (B24) 

r— ^0 



This justifies the counting rule of the number of regular 
solutions at the origin in Sec. IIVI 
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